Pauli-Limited Superconductivity with Classical Magnetic Fluctuations 
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We examine the effect of classical magnetic fluctuations on the phase diagram of paramagnetically- 
limited two-dimensional superconductors under a Zeeman magnetic field. We derive the free energy 
expansion in powers of the superconducting order parameter and analyze the character of the normal- 
superconducting transition. While the transition is of the second order for all temperatures in the 
absence of magnetic fluctuations, we find that proximity to magnetism drives both the transition 
into the uniform state and that into the modulated (Fulde-Ferrell-Larkin-Ovchinnikov, FFLO) state 
to first order at intermediate temperatures. We compute the thermodynamic signatures of the 
normal-superconducting transition along the upper critical field. 
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I. INTRODUCTION 

Since the discovery of the Meissner effect, there has 
been ongoing study into the interplay of superconductiv- 
ity with magnetic fields. In type II superconductors, su- 
percurrents, which arise from the coupling of the Cooper 
pair momenta to the vector potential of the magnetic 
field, cause the system to return to the normal state at 
the orbital-limited upper critical field Hc2,orb- If, how- 
ever, the dominant coupling of spin-singlet superconduct- 
ing electrons to the field is via the Zeeman effect, the 
transition occurs at the Pauli-limited field Hp where the 
condensation energy equals the gain in energy due to spin 
polarization of the two electrons in a Cooper pair. In ad- 
dition, in pure systems in the vicinity of this field, the 
superconducting state can exhibit spatial-modulation of 
the type predicted by Fulde and Ferrell, and Larkin and 
Ovchinnikov (FFLO) at low temperature and high mag- 
netic fields^'^. 

No clear examples of FFLO superconductivity have 
been found; however, early experiments^""— on the 
layered^, heavy-fermion superconductor CeCoIus, tenta- 
tively identified the low-temperature, high-field (LTHF) 
superconducting phase as a possible realization of the 
FFLO state. The heavy mass and large value of the ra- 
tio between the estimated orbital critical field and the 
Pauli-limiting field y/2Hc2,orb/ Hp ~ 3.5,— suggest strong 
paramagnetic limiting. It was also proposed that the un- 
usual field dependence of the vortex lattice form factor 
of CeCoIus is due to Pauli-limiting effects^ii^. 

Some experimental features of the transition into the 
LTHF phase are not fit by the established theories of 
Pauli-limiting. For example, it has been established 
experimentallj*^ that the transition from the normal (N) 
to superconducting (SC) state in CeCoIus is first order 
into low-field state (T < Tfflo), and remains first order 
at low temperatures up to Tq > Tpp^o- This is in sharp 
contrast to the conventional theory in two dimensions 
(2D) that finds a second order N-SC transition along the 
entire critical field line Bc{T)p^^, and motivates our cur- 
rent study. 

Under purely Zeeman field, B, when the electron 



spins couple to the field but the orbital coupling to 
the vector potential is irrelevant, at low tempera- 
tures the N-SC transition is into an inhomogeneous 
state^'^. In most cases the amplitude modulated state, 
A(r) — Aq s\Ti{Qx) (LO state), is favored compared to 
the purely phase-modulated FF state, A(r) — Aq e''^^. 
In the conventional analysis for isotropic, s-wave, super- 
conductors the FFLO transition is 2nd order in 2D and 
1st order in S D^^i^^" — . In superconductors with nodes, 
such as d-wave CeCoIug^ii^""— , the FFLO transition is 
2nd order in both 2D^^ and 3Di^. The transition to a 
combined vortex and LO state is also expected to be sec- 
ond order—. 

Theoretically, under several conditions the transition 
from the normal to the SC state may become first or- 
der. For purely Zeeman coupling this can happen a) due 
to strong Fermi-liquid enhancement of the magnetic sus- 
ceptibilit y^^'^^ , or b) due to impurity scattering in the 
resonant limit^^. In the presence of both paramagnetic 
and orbital effects it was argued in Refs. [22.23] that the 
transition may also become first order in an intermediate 
temperature range. 

In this paper we show that in systems with enhanced 
magnetic susceptibility (as compared to the dimension- 
less Pauli susceptibility for typical metals xp ~ 10~^), 
such as some heavy fermion materials, first order N-SC 
transitions under purely Zeeman magnetic field may nat- 
urally emerge over a part of the phase diagram. Impor- 
tantly, the transition is first order both for the uniform 
and for the FFLO state over part of the temperature 
range. We consider the critical field Bc{T) and discuss 
how the magnetic fluctuations affect the order and ther- 
modynamics of the transition. 

The rationale for inclusion of such fluctuations is as 
follows. Experiments convincingly show that CeCoIus 
is in proximity to a magnetically ordered state^-i^""— . 
In CeCoIus the /-electron spins are not fully Kondo 
screened by the onset of the superconducting orderi2^, 
and the entropy of the remaining spin fluctuations is re- 
leased at the superconducting transition'^'^. The specific 
heat jump upon entering the superconducting state at Tc 
in zero field is AC/jT^ ~ 4.5 (where 7 is the Sommerfeld 
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coefEent), more than three times the s-wave BCS value 
of 1.4S2.. Guided by this observation, Kos, Martin, and 
Varma^ considered a Ginzburg-Landau model of compe- 
tition between thermal (Gaussian) magnetic fluctuations 
and superconductivity, and were able to explain the large 
specific heat jump at Tc{B — 0). 

We employ similar methods, with more microscopic 
considerations, to include an applied magnetic field. To 
simplify our analysis, we consider only the Zeeman cou- 
pling since it is largely responsible for the salient features 
in the phase diagram {e.g., the first order transition). We 
explain our results in context of experiment for CeCoIns 
and emphasize that our approach is generally applicable 
to Pauli-limited systems with thermal magnetic fluctua- 
tions. 

The rest of the paper is organized as follows. In Sec. II 
we review the physics of superconductivity in the para- 
magnetic limit and our methodology of determining the 
transition line. In Sec. Ill we extend the formalism to ac- 
count for the magnetic fluctuations, and use this method 
to obtain the results presented and discussed in Sec. IV. 
We conclude by placing our results in the context of ex- 
periment and theory on Pauli-limited superconducting 
systems. 



II. SUPERCONDUCTIVITY IN THE 
PARAMAGNETIC LIMIT 

A. Model and approach 

In the paramagnetic limit the second order transition 
into the uniform superconducting state, A(r) = Aq, 
becomes first order below a characteristic temperature 
Tp ~ 0.56Tco where T^q is the transition temperature in 
zero field^i. At T = superconductivity is destroyed 
when the energy of the polarized normal state equals the 
superconducting condensation energy (the Clogston'^-- 
Chandrasekhar'^^ limit). This occurs at the Pauli field 
Hp = A/{\/2^), where — 3/15/2 is the electron 
magnetic moment, /.ip is the Bohr magneton, and g is 
the conduction electron g-factor. The microscopic pair- 
breaking occurs as the Zeeman field increases the en- 
ergy of the spin-singlet with respect to spin-polarized 
s — 1. An alternative to this uniform superconductiv- 
ity is the pairing of the electrons with opposite spins and 
the same energies, which now have momenta differing by 
Q ~ fj,B /(hvp), where vp is the Fermi velocity. The fi- 
nite center of mass momentmn of the Cooper pairs leads 
to a spatial modulation of the order parameter^'^ and 
allows superconductivity to survive at fields above the 
Clogston-Chandrasekhar limit. 

The exact structure of the modulated state in s- 
wave systems is still not well established. Generally the 
amplitude-modulated LO state is lower in energy than 
the current-carrying FF statei^. In the absence of spin- 
orbit coupling, the direction of Q in real space can be cho- 
sen arbitrarily, and a superposition of plane wave mod- 



ulations along different directions may yield yet lower 
energji^Siil. In systems with unconventional, such as d- 
wavei^^iSi^i^, gap symmetry the modulation is prefer- 
entially along either the nodal or anti-nodal orientation, 
depending on both the temperature and the purity of 
the samplei^i2i, with the LO state always more advanta- 
geous. 

Analysis of the FFLO states is often carried out within 
the Ginzburg-Landau theory, expanding the free energy 
of the superconducting state in both the amplitude and 
the gradient of the order parameter—. Such an expan- 
sion is justified in the immediate vicinity of Tp, but its 
region of validity is very narrow. The modulation wave 
vector increases rapidly along Tc{B) < Tp and becomes 
comparable to the inverse of the SC coherence length, 
CcT^ = [Tivp /2TrTco] ^, rendering the gradient expansion 
invalid. Previously we reported the results of a brief anal- 
ysis of the N-SC transition for a superconductor with 
magnetic fluctuations under Zeeman field using such a 
gradient expansion—. The results were suggestive of the 
onset of the first order transition near Tp. The limita- 
tions of the gradient expansion prevented us from reach- 
ing detailed conclusions, and motivated our present work. 
Below we expand the free energy only in powers of the 
order parameter and retain the full wave-vector depen- 
dence of the expansion coefficients, thus removing the 
deficiencies of the gradient expansion and presenting a 
theory valid down to low temperatures. This allows us 
to analyze the details of the phase diagram not accessible 
with the gradient expansion. 

We begin with the mean-field Hamiltonian 

(1) 

q.k 

where cr = ± denotes the orientation of the electron with 
the spin along/opposite to the field direction, ek is the 
band energy measured with respect to the chemical po- 
tential, and ek.cr = Ck + ufiB. In Eq. ([1]), |A| is the 
strength of the pairing interaction, 3^(k) is a normalized 
basis function that transforms according to an irreducible 
representation of the crystal point group and describes 
the gap symmetry, and k denotes position on the Fermi 
surface (FS). 

We assume, for simplicity, a separable pairing interac- 
tion, so that the spin-singlet order parameter is V'(k, q) ~ 
3^(k)Aq, with the amplitude Aq self-consistently deter- 
mined from 

Aq = -|A|E3^(k)(ck+q,+c-k.-) (2) 

k 

where (■ ■ • ) indicates thermal average. Uniform super- 
conducting states have the single non-vanishing Fourier 
component with q = 0, while modulated states corre- 
spond to one or more components with q 7^ 0. Since 
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Eq.(l2]) has to minimize the free energy, it determines, at 
the mean field level, the Landau expansion of the free 
energy density Fl in powers of Aq, 

{q.} 

+ 51 T'qi'--q4^qi^q2^q3'^q4^qi+q3,q2+q4 

{q.} (3) 

+ E ^qi. -qe^qi^q^^qsA^.Aq^A*^ 
{q.} 

^ '^qi+q3+q5,q2+q4+q6 • 

The coefficients of this expansion are combinations of 
the normal state Green's functions as described in Ap- 
pendix The summation over {qi} in Eq. ([3]) includes 
all possible combinations of the allowed Fourier compo- 
nents of A(r): qi = for a uniform gap amplitude, single 
mode qi = Q for the FF modulation, and q^ G {Q, — Q} 
for the LO phase. In the following we restrict ourselves to 
the comparison of the free energies of these three phases, 
finding the one most energetically favorable and the cor- 
responding wave vector Q. 

Based on the observation of the quasi two-dimensional 
Fermi surface in the 115 family^"—, we use a model of a 
2D circular Fermi surface. We use the azimuthal angle, 
0, to parameterize the position on the Fermi surface, and 
we choose 3^(k) = y{0) = 1 and y{e) = ^2 cos 261 for s- 
and d-wave gaps, respectively. 

We determine the phase transition line Bc{T) by 
finding, at a given temperature, T, the highest Be of 
the three phases we compare. In each phase we find 
Bc{T) = max(_Bc(T', q)) by unrestricted maximization 
with respect to the modulation wave vector. We intro- 
duce the dimensionless energy density, / = Fl/NpT"^^, 
where Np is the 2D normal state density of states at 
the Fermi level. We also introduce the dimensionless 
amplitude 5q = Aq/Tcq where Aq is the SC gap am- 
plitude and Tco is the mean-field transition temperature 
at -B = in the absence of magnetic fiuctuations. We 
set fcs = = 1 throughout the paper. The reduced tem- 
perature and magnetic field are given hy t = T/Tcq and 
b = fiB / {2TrTco) respectively. 



B. Uniform superconducting state 

For the uniform state, Aq = Ao(5q.o, we find 

/„(r, B) = a„|5o|' + lu\5o\^ + i^u\Sof , 
with the coefficients determined from Eqs. 



In (t) + Re 



^ { - +1- 
2 t 



8 (27rt)2 

1 {\ym')ps 



Re 



192 (27rt) 



(4) 

, (5a) 
(5b) 
(5c) 



Here ^ (v]/(")) is the digamma (nth order polygamma) 
function, and {■■■)fs — J d9/{27T). For the s- and d- 
wave symmetries of the gag, Eq. ()A12p . our coefficients 
agree with those in Refs. [isHTgj . 



C. FF state 

For the spatially-inhomogeneous superconducting 
state, the coefficients in Eq. ([3]), depend on the direc- 
tion of modulation. Since the modulation wave vector 
Q ^ CfT^ ^ kp, for two particles at locations 9 and 
t: + 9 on the Fermi surface, there is an energy mismatch 
vf ■ Q — vfQ cos{9 — 9q), where 9q is the modulation di- 
rection with respect to the crystalline a axis. This energy 
mismatch enters in Eq. (jA8l) with q^ = Q. 

Recall that the polygamma functions in Eq. ([5]) orig- 
inate from the summation over Matsubara frequencies, 
and that their argument is determined by the energy 
mismatch of the particles in the Green's functions in 
Eqs. (IA8P - (jA10[) . Consequently, the coefficients of the 
free energy expansion in the FF state are given by 
the same polygamma functions as for the uniform case, 
Eq. ([5]), but with the arguments reflecting the energy dif- 
ference ^B + vp -Q. Hence in the expansion fpp(T, B) — 
app\Sa\'^ + 7_ff|'5o|'* + i^ff\So\'^, we find 



apF = In (i) - * ( ^ 



Re(^\yi9)\H(^\+^^^ 



(6) 



FS 



where q — qcos{9 — 9q) and q — ^o(9/2. Similarly, ^pp 
and vpp are given by expressions identical to Eqs. (jSbp 
and (|5cp under the replacement b ^ b + q and averaging 
both the digamma functions and the basis functions y{9) 
together over the Fermi surface. 

It follows that for any anisotropic superconductor the 
direction of the modulation and the shape of the gap 
cannot be separated. For a two-dimensional d-wave su- 
perconductor that we consider, the modulation along the 
nodal/antinodal direction is preferred in a pure material 
above/below T ~ O.OeTccr^'^'^, although as the impu- 
rity scattering is increased modulation along a node be- 
comes favorable even for T < O.OGTco^^: Therefore, below 
we focus on the modulation along the gap nodes. 



D. LO state 

For the Larkin-Ovchinnikov (LO) state, the quadratic 
component in Eq. Q includes two terms identical to 
Eq. ([H but summed over q — ±Q with A±q = Ao/2. 
Both terms for LO are identical when averaged over the 
Fermi surface, hence ulo = app/2. Thus the second or- 
der transition line, Be, determined from a = 0, is identi- 
cal for both the FF and LO phases. The relative stability 
of the FF and LO phases is determined by comparing the 
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quartic coefBcients 7f_f and 7lo at the transition, with 
the smaller of the two corresponding to the thermody- 
namically stable SC state because fsc — In = — Q!^/(27)- 
The quartic coefficient, 7lOj is obtained by summing 
the six terms in Eq. (IA9p with € {Q, — Q}, subject to 



the constraint 5a 



This yields 



ILO 



128 



(7) 



FS 



where ujn,b = tin + i) + ife. Twenty distinct terms con- 
tribute to the sixth order Landau coefficient which be- 



comes 



\n=0 



2048^4 ^2 ,) (9g-2 + ^2 ,) 



FS 



We can obtain the gradient expansion of the free en- 
ergy by expanding Eqs. ([B])-® in powers of q with the 
corresponding y{9). The resulting Ginzburg-Landau ex- 
pansion coefficients are identical to those obtained for s- 
and d-wave SC in Refs. [Tsl and [l^, respectively. Below, 
however, we retain the full g-dependence of the expansion 
coefficients to examine the transition line at low-T where 
Q ~ ^(^^ and the gradient expansion s ^^i^^ fail. 



E. Determination of Bc{T) 

For each phase, with the free energy density written as 

/ - a{t, b, q)\6o\^ + 7(i, b, q)\So\^ + v{t, b, g)|5o|' , (9) 

we determine the critical field Bc(T) and the optimal 
modulating wave vector qq. We allow for possible second 
and first order transitions, and compare the results to 
determine the order of the physical transition. 

The second order transition field at fixed t is the max- 
imal value of be (with respect to q) for which a{t, bc,q) = 
0, and j(t,bc,qo) > 0. The corresponding optimal qo de- 
termines whether the transition is into a uniform (go = 0) 
or modulated state. In the vicinity of the transition line 



a{t,bc,qo) a'{tc-t) 



2j{t,bc,qo) 2j{tc,bc,qo) ' 



(10) 



where a' = da{t,bc,qo)/dt\t=t^- With this value we can 
compute the free energy difference between the normal 
and the superconducting states and therefore determine 
the thermodynamic properties such as the specific heat 
jump at the transition, see below. 

In the region where 7 < 0, the first order transition 
occurs once the minimum in the free energy shifts dis- 
continuously to 7^ before the quadratic coefficient a 
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FIG. 1: (Color online.) Optimal wave vector go (fee (r)) (solid) 
and quartic coefficient ^{Bc{T)) (dot-dashed) obtained by un- 
restricted maximization of bc{t) for LO state. Quartic co- 
efficient for small-g expansion (dashed) is also shown. Up- 
per/lower panels are for s and d-wave (g along node) gaps. 



changes sign^*^. This happens along the line defined by 

j\t, 6„ q) - 4a(i, 6„ qHt, 6„ g) = , (11) 
where the new minimum first appears at 



\X |2 "f 



(12) 



We locate the first order transition at a given t by 
unrestricted maximization, with respect to q, of the field 
be that satisfies Eq. ((TT|) . 

At each temperature, we locate the maximal field 
bc{q) = bc{qo) for which the coefficient a{t,bc,qo) — 0. 
If we find "f{t,bc,qa) > 0; the transition is second order. 
If "f{t,bc,qa) < 0; we maximize be for Eq. (jlip . checking 
that 7 remains negative and that the free energy remains 
bounded from below (z^(i, 6c, <Zo) > 0). The first and sec- 
ond order transition lines meet at a critical point t* where 
j{t*, be, qo) = 0. For cZ-wave gap, we compare the critical 
fields for nodal and antinodal orientations of the modu- 
lation wave vector Q and verify that modulation along 
gap nodes is preferred above T ~ O.OGTcp. 

As shown in Fig. [T] for LO modulation, the quartic 
coefficient remains positive and the transition is second 
order on both sides of tp. For comparison, we also plot 
the results obtained from the small q = expansion 

f^{ao + a2q^ + aiq^)\5o\^ + (70 + l2q^)\5o\'^ + va\5o\'^ 

of Eq. (ini), where the coefficients are found by expand- 
ing Eqs. (IA8l ) and (jA9p . and only even powers of q ap- 
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pear since the system is isotropic. Since each subse- 
quent term in the expansion contributes an extra qG^^ 
in the Matsubara summation of Eqs. (jA8| )- ()A9|) . we have 



a„ cx 7„_2 oc Re(^'(")(i + ij)) for n > 2. Conse- 
quently a2 and 70 change sign at exactly the same tem- 
perature, tp = Tp/TcQ, and the modulated state with 
Qq = — Q!2/(2a4) emerges at lower T via a second order 
transition described by 



/ = "0 



A.] 

4a4 / 



70 



2a4 / 



(13) 

with the renormalized quartic term positive. The results 
obtained from an examination of Eq. (|13p are identi- 
cal to those discussed in Refs. ^M^- As is seen from 
Fig. [1] below tp the transition is into modulated state 
with the wave vector that reaches values of qo ~ 0.2^0^^ 
and higher. The failure of the gradient expansion is man- 
ifested in the significant discrepancy between the val- 
ues for the quartic coefficient at the optimal wave vec- 
tor within the gradient expansion and in the full evalua- 
tion, shown in Fig. [T] Hereafter, we rely on unrestricted 
maximization of bc{t) with respect to q to determine the 
critical field. 



III. MAGNETIC FLUCTUATIONS 
A. Magnetic fluctuations 

If we are close to a magnetically ordered state, soft 
magnetic modes exist in the system. In the continuum 
limit the fluctuations of the magnetization field M(r) are 
described by the Gaussian free energy 



[M] 



1 J drx''MirY 



(14) 



Here we have ignored the the momentum dependence of 
x(q), assuming the that the momenta relevant for super- 
conductivity are Q ^ n/a (where a is the lattice 
spacing) . We do not discuss here the role of the (possibly 
singular) antiferromagnetic (AFM) fluctuations in medi- 
ating the superconducting pairing: this role can only be 
addressed within the framework of specific microscopic 
theories^"—. Our task is to consider the competition of 
superconductivity and the long-wavelength fluctuations 
of the magnetization, whether uniform or staggered. Al- 
though the susceptibility is enhanced (x ^ Xp)i the sys- 
tem is not close to ferromagnetic order (x <^ 1), hence 
we do not distinguish between B and the applied mag- 
netic field H for the rest of the paper. In the same spirit, 
we ignore M • H in the magnetic free energy since its 
contribution to the averaged free energy is a factor of x 
smaller than the corrections we consider. 

In Ref. [3^ the susceptibility was taken to be 
temperature-dependent, in agreement with experi- 
ment?^, x{T) = XoTsf/{T + T,f) where T^f is a charac- 
teristic energy scale for the low-energy spin fiuctuations. 



While we make use of this expression to make contact 
with Ref. [3^ , our main results are qualitatively the same 
for a temperature independent x of the same magnitude. 
We also ignore the field dependence of x- Finally, we 
do not account for the quantum fluctuations of M and 
consider only thermal fluctuations of the magnetization. 

Recall that our goal is to investigate the effects of long- 
wavelength magnetic fluctuations on the N-SC transition. 
We include the competition between magnetism and su- 
perconducting orders via the lowest order term allowed 
by symmetry in the free energy expansion. 



J-,,,m[A,M] = | j dr|A(r)|2M(r)2, 



(15) 



where the coupling constant 77 > makes coexistence of 
the two orders unfavorable. In a simple system it would 
be possible to determine 77 from microscopies, by expand- 
ing B = H -I- 47rM in each Green's function in the pow- 
ers of the fluctuating magnetization and introducing the 
correlator (M(r)M(r')) that is prop ortional to the sus- 



ceptibility, in analogy with Ref. [5J|. Such an expansion 



produces an M |A| term. In a complex system with 
/-electrons we cannot determine the coefficient of this 
term microscopically, and we use Eq. p5|) with a phe- 
nomenological parameter rj to explore the salient features 
of the model. Our consistency checks on the choice of 77 
are the magnitudes of the jump in Aq and of AM/M 
across the first order transition line. We find maximal 
Ao(Tc) < 0.3A(0) and AM/M « 1-5% everywhere 
along the first order transition line. These values are 
moderate, hence our choice of rj is physically reasonable. 

To verify the ubiquity of our results, we examined the 
coupling of the SC order parameter to higher order terms 
in M(r)2 and its gradients, e.g. |A(r)p|VM(r)p, and, 
within the small-g approximation, to gradients of the or- 
der parameter itself, e.g. |VA(r) • M(r)p. We checked 
that, while these various couplings renormalize the tran- 
sition temperature, they do not introduce new features 
into the phase diagram. 

To derive the effective theory for the superconducting 
order we integrate out the magnetic fluctuations from the 
partition function 

Z = exp [-{Tl + J'sc.M + ^m)IT] = e-^^'^^sc,M 
where Fl = / dPvFL. We obtain the total free energy 



where the magnetic contribution is 



■M 



2?[M(r)] exp 



M, 



(16) 



(17) 



and I?[M(r)] indicates integration over all possible con- 
flgurations of magnetization. The integral is Gaussian 
in M, hence we compute it analytically and expand in 
powers of |Ap to obtain the corrections due to magnetic 
fluctuations to the expansion coefficients in Fj^. Below 
we address these corrections in each of the three phases 
we consider: uniform, FF, and LO. 
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B. Uniform superconducting state 

Integrating out the fluctuating magnetization for a uni- 
form order parameter, Aq, is straightforward. We work 
with the Fourier components of the magnetization, Mk, 
and restrict the sum 



IMkP. 



(18) 



k|<fcc 



to one-half of fc-space since Mk — M* for real M(r). 
Therefore from Eq. ((T7l) we have, after Gaussian integra- 
tion over both real and imaginary parts of Mk, 



■M = 



n 

k|<fcc 



l + 77X|Ao 



(19) 



where now the product is taken over all k up to the cutoff 
of the order of the lattice spacing |kc| = tt/1, and d is the 
dimensionality of magnetization vector M. 

Neutron scattering^ measurement of the dynamic spin 
susceptibility in CeCoIns shows evidence of spin fluc- 
tuations, and light Cd-doping^^ induces AFM order at 
Qafm = (-5, .5, .5)22.. Sister compound CeRhlus ex- 
hibits AFM order at Qafm = (-5, .5, .297)^^ which is 
stable under pressur o^^'^^ , before SC preempts AFM or- 
der at P 2GPa^^. Furthermore, the pressure depen- 
dence of Tc and T/v for CeRhlus and Cd-doped CeCoIus 
is nearly identical^^ suggesting that the SC and magnetic 
orders in both are closely related. Hence we conclude that 
CeCoIns is in proximity to 3D magnetic ordering, and we 
take d = 3 for the purposes of this paper. 

The corresponding contribution to the free energy is 
J-'Af(Ao) — —TlnZsc,M- Subtracting the average mag- 
netic contribution to the normal state energy, Tm{^o = 
0) we find an additive contribution to the superconduct- 
ing free energy density 



uni.M 



(0) 



= |td E l^^(l + ^^l^on (20) 



2LD 

= ^^ln(l + r;x|AoP) 

where is the volume. The last line of Eq. (f20| follows 
from '^f.a ~ a{L/l)^ where a does not depend on k. 
Expanding this contribution in powers of |Ao| for our 
2D superconductor [D — 2), we find the renormalized 
coefficients of / = F/{NfT^q), Eq. ®, 

3 . 



Ctu + -zt 



2 Npf 



lu = In 



3^_T3 



4 Npl"^ 
2 NfP ' ^ ^ ' 



(21a) 
(21b) 
(21c) 



where a„,7„, and fu are given in Eqs. (l5al) - ([5c)) . 

Since in the FF state, A{x) — Aqc^'^^, only the phase 
of the order parameter is modulated, the coupling be- 
tween the magnetization and the superconducting order, 
Eq. (jisp . has exactly the same form as in the uniform 
state. Hence the renormalized expansion coefficients are 
obtained from Eqs. (j21[) by a direct substitution of apF, 
jFF-, and lypF for Uu, 7u, and Vu-, respectively. 



C. Modulated LO state 

In the LO state, in addition to the order parameter 
A{x) = Aq cos{Qx) competing with the average magne- 
tization, the amplitude modulation couples the magnetic 
fluctuations at wave vectors differing by 2Q. Therefore, 
the magnetic contribution to the free energy is 



J'scAI + J'M — E 
|k|<fce 



1 



1 



X 



Vl^ol' |Mk 



--77|AopMk-(M*+2Q + M^_2Q) 



(22) 



After integrating out the fluctuations, the contribution 
to the superconducting free energy density relative to the 
normal state becomes 



i'LO.M - 2^ 



1 



ln( l + -ryx|AoP 



(23) 



ln(l--,7V|Ao|^ + -r7V|Aor 



The first term differs from its counterpart in Eq. ((20)) 
by the factor of 1/2, arising from the spatial average of 
cos^{Qx). The second term arises from the mode- mode 
coupling terms in Eg. (1221) and is derived in Appendix IB] 
Under expansion in Aq, it only contributes to the fourth 
and sixth order terms in the free energy, and we obtain 



aLO = aLO + ^*]^^'?^(^)' 

— -'-cO 2 2/rTi\ 



(24a) 
(24b) 
(24c) 



where aLOjlLO, and vlo are given in Sec. Ill Dl 

Comparing Eqs. (|21l) and ([24l) . we see that the free 
energy expansion depends on t] and x only through their 
product r]x- Thus, for subsequent analysis we define a 
dimensionless coupling parameter 



3 Tco 3 

V — ttV — —IrCl-L FV 



(25) 



where the characteristic temperature Tp = {NpP) ^ is 
of the order of the Fermi temperature in the system. We 
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also define a dimensionless parameter based on the ex- 
perimental fit of x{T) = XoTsf/{T + Tsf) 

Tsf 



X = Xo 



TcQ 



(26) 



With these parameters, the renormalized quadratic coef- 
ficients in Eqs. (|21al) and (I24al) become simpler, e.g., 

t , , 

a« = a« + fix——— ■ (27) 

t + tsf 

The renormalization of all other Landau coefficients is 
determined by the product fjx, and, in simplifying the 
fourth and sixth order terms in Eqs. ([2T|) and (f24| . we 
introduce the parameter tp — Tp/TcQ. 

We note that the dimensionality of the magnetization 
vector M enters the Landau coefficients as a prefactor 
of the coupling parameter 77. Throughout this paper, we 
take d — 3. Using a different value for d simply decreases 
the magnetic fluctuation contributions in Eqs. (j21[) and 
(l24t by a factor of d/3. For example, taking d — 2 only 
requires that we use 877/2 to obtain the same results {e.g., 
Tc) as for rj and d — 3. Hence, we proceed with our choice 
d — 3 without any loss of generality. 



D. Choice of energy scales and parameters 

The exchange of entropy between the magnetic fluctu- 
ations and superconductivity reduce the zero-fleld tran- 
sition temperature from the unrenormalized Tcq to the 
experimentally observed Tc{ri > 0) as determined from 
the instability condition 

= au = In (^^^ + It.Tpijx ■ (28) 

The extra entropy is released in a speciflc heat jump that 
exceeds the BCS value. 



QT2 



[a'iv)] 



27(77) 



(29) 



Taiv) 



where /(?/) is the dimensionless free energy for the given 
coupling, 77, and a' — da{T)/dT. Without magnetic 
fluctuations, BCS mean field theory predicts for s-wave 
gap AC/Cn — 12/7^(3) sa 1.43 and for d-wave gap 
AC/Cn = 8/7C(3) « 0.95 at T^q. Here Cn is the nor- 
mal state specific heat, and ^(3) w 1.202 is the Riemann 
zeta function. Measuring the jump relative to the s-wave 
value, we find for B = 



1.43Cn/T, 



'1 



lT,Tp7^{x + T,x')y 



cO 



mow 



7C(3) ^c^F'l X 



(30) 



where = 7;(77) and x' = dx{T)/dT. Using Eq. ^ 
to eliminate rj we find 

l + ^^-^ln(^))' 

(31) 



AC/Tc 



1.43C./T.0 (|y(e)n-gf§ln^(^)' 



in zero field. We discuss the field dependence of AC/Tc 
in Section HVBH 

From the experimentally measured behavior of the sus- 
ceptibility, specific heat jump AC/Tc{B = 0), and T^. one 
can estimate T^q provided a reasonable guess about the 
value of Tp can be made. For our purposes, we take 
Tc = 2.3K, Tp = AOK (the Kondo coherence tempera- 
ture for CeCoIng^), and the dimensionless xo ~ lO""* 
(presented in units of emu/g in Ref. [1^). We follow the 
example of Ref. [s^ and set T^f = 1.5K. With this choice 
Tsf < Tc, and we examine the effects of x which varies 
substantially with temperature below Tc. Experiment, 
however, suggests a weaker temperature dependence of 
X(r) with Tsf « 3.5T^. Therefore, we verify that our 
general results are independent of the details of x by com- 
paring this case with the analysis for constant susceptibil- 
ity. For our chosen energy scales, we solve Eq. (I5T]) with 
AC/Tc = 3AC/Tco. This gives for s-wave Tco = 
and fjx — 1.6 and for d-wave Tco = 9.27 K and fjx — 2.3. 



IV. DISCUSSION AND RESULTS 

Using the formalism outlined above, we are now in the 
position to investigate the changes appearing in the tran- 
sition lines of the superconductor coupled to the magnetic 
fiuctuations. In the following we set Tp, Tsf, and TcO as 
described at the end of the previous section. We adjust 
the coupling 77 to the magnetic fiuctuations as well as the 
temperature dependence of the magnetic susceptibility. 
We first address the nature of the transition along the 
Bc{T) line, and then consider the thermodynamic signa- 
tures of these transitions. 



A. Normal to superconducting transition in a 
magnetic field. 

Quite generally coupling to magnetic fiuctuations sup- 
presses the transition temperature, since, as is clear from 
Eq. ((T5)) . the finite thermal average of M^(r) makes 
the appearance of superconductivity energetically costly. 
This is also evident from Eas. (|24ap and (|21al) . which 
show positive additive contribution to the quadratic co- 
efficients in the Landau expansion. In the absence of the 
field, when a{T) = -InT/TcO, it follows from Eq.(| 
that the transition temperature Tc satisfies 



Tc_ 

TcQ 



= e-i^'^-^^(T^'> = cxp 



B=0 



'Tc + T 



(32) 



where in the last step we explicitly invoked the temper- 
ature dependence of the susceptibility. For small 77% the 
linearized form of this equation coincides with that used 
in Ref. [361]. 

At the same time the results for the quartic coeffi- 
cient, Eas. (j24bl) and (j21b|) show that it is renormalized 
downward by the magnetic fluctuations. Since the sign 
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0.0 0.25 0.5 0.75 1.0 

TIToo 



FIG. 2: (Color online.) The normal-to-superconducting tran- 
sition in s-wave (upper) and d-wave (lower) superconductors 
under a Zeeman field. Modulation is along gap node for d- 
wave. Magnetic fluctuations, t), modify the Bc2 transition, 
showing 4 distinct regions (for increasing T): 2nd order into 
modulated, 1st order into modulated, 1st order into uniform 
and 2nd order into uniform states. 



of this term controls whether the transition is of the sec- 
ond or first order, it seems possible that the order of the 
transition may change as the strength of the magnetic 
fluctuations increases. 

Figures [2] and [3] show that this is indeed the case: cou- 
pling to magnetic fluctuations opens a region of first or- 
der transition from the normal to both uniform and the 
modulated superconducting state. This finding is a ma- 
jor conclusion of our work, and qualitatively fits with 
the behavior of CeCoIn5 where the transition becomes 
first order below Ti « IK—, while the putative FFLO- 
like phase does not occur until a lower temperature^. 

To understand this behavior recall that in the absence 
of fluctuationsi^iii^i^ the quartic term of the Ginzburg- 
Landau expansion for the uniform superconducting phase 
changes sign, ^u{Tp) = exactly at the point along the 



0.3 I . ^ . ^ r 




T/T^ 

cO 

FIG. 3: (Color online.) Transition lines for s- and d-wave 
with nodally-oriented Q. Upper curves are in the absence of 
fluctuations, and lower curves are for fjx = 0.5. 

Bc{T) line (at temperature Tp) where the modulated 
phase, reached via a second order transition, aLo{Tp) = 
0, becomes allowed. Coupling to the fluctuations in- 
creases oiLO a-iid lowers 7„ ensuring that the first order 
transition in the uniform state occurs at higher temper- 
ature than that where the modulated phase can form. 

As is seen from Figs. [2] and [3] the region of the first or- 
der transitions widens as the fluctuations become softer 
{x increases) or compete more strongly (ry increases) with 
superconductivity. There we define the temperatures t*p 
and t^Q where the second order transitions into the uni- 
form and LO modulated superconducting states respec- 
tively become first order. 

Since we use the expansion in powers of 5q we can 
only estimate the location of the first order transition 
line away from the critical points at which the transition 
becomes second order. However, since the jump in 5q 
across the first order transition is modest {e.g., for s-wave 
gap, 6o{tc) < O.3(5o(0) with (5o(0) = ttc-t^ « 1.76) this 
estimate is quite reliable. We denote by Ilo our estimate 
of the temperature along the bc{t) where the first order 
transition lines into the uniform and the LO phases meet. 
For t < Ilo the transition (first or second order) is into 
the amplitude-modulated phase, while for t > tj^o it is 
into a uniform phase. In the absence of fluctuations, of 
course, t*p = t^Q = Ilo = tp- 

We find that for s-wave order the region of the first 
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cO 

FIG. 4: (Color online.) Quartic Landau coefficient evaluated 
along Be (T) for s-wave and d-wave symmetries with LO mod- 
ulation for t < tp and uniform state for t > tp. 



order transition, for the same values of the coupling and 
magnetic susceptibility, is wider. This can be qualita- 
tively explained by examining the quartic Landau coefH- 
cient for both symmetries in the absence of fluctuations, 
shown in Fig. H) In the vicinity of Tp, the coefficient 
7 is numerically smaller for an s-wave order parameter 
than for d-wave, both on the uniform and the modulated 
(with the wave vector Q yielding maximal for each 
symmetry) side of the transition. Hence it is easier to 
drive an s-wave system to first order transition. 

Note that for d-wave SC we find that the modula- 
tion of the order parameter along the gap nodes is sta- 
bilized even below T — O.OGTcOj where, in the absence of 
fluctuations, the anti-nodal direction would be more ad- 
vantageousi^^"^. The anti-nodal modulation still gives 
a lower free energy at very low temperatures, below a 
threshold that depends on the parameter f)x, but that 
occurs far from the first order transition range that is 
our focus here, and therefore for the rest of this paper, 
we discuss only d-wave SC where Q is oriented along a 
gap nodes. 

The key finding of the region of the first order tran- 
sition does not depend on the exact temperature depen- 
dence of x(^)- For comparison, we also considered the 
constant susceptibility xi = x(2^c) so that, for a given 
coupling strength fj, we obtain the same T^- In Fig.O we 
compare the critical field and order of transition for x{T) 
and constant xi- Since fjxi < fix{T) for all T < Tc, su- 
perconductivity is suppressed less and Bc{T) is higher for 
constant susceptibility. However, in both cases the prod- 
uct f)x < 1, and the magnetic fluctuations have a larger 
effect for constant susceptibility than for x{T) on the 
fourth Landau coefficient where fjx enters quadratically. 
Thus, the N-SC transition is first order over a wider tem- 
perature range for constant susceptibility. While both 
the exact temperature range of first order transition and 
the degree of Bc{T) suppression depends on the temper- 
ature dependence of %, the presence of these effects is 
independent of the details of the susceptibility. Further- 




FIG. 5: (Color online.) S'-wave bc{t) for both constant x(i) 
and Xi = x(^c) for tc = 0.76. Upper and lower curves are for 
Xi and x(*)i respectively. The critical field is suppressed less 
for xi since fjxi < VXit) at t < tc. As f/x < ^ for both cases, 
the region of first order transitions is larger for xi- 



more, the thermodynamics of the transition are similar 
for both xi and xC^) where the only significant differ- 
ence is the low-T behavior of the specific heat jump for 
d-wave as discussed below. 

We note that our results agree with those obtained 
from a small-q expansion of the free energy functional^. 
In Fig. ini we compare the results obtained from the fully 
q-dependent Landau functional and the small-g approx- 
imation by plotting the optimal wave vectors found via 
each method. Each model predicts first and second or- 
der transitions into both the uniform and modulated SC 
states. Hence, our current model supports our prelim- 
inary results^ while allowing us to examine the upper 
critical field beyond the limitations of a small-g approxi- 
mation. 

Our main conclusion so far is therefore that coupling 
to thermal magnetic fluctuations drives the transition 
from the normal to superconducting state first order in 
the vicinity of the onset of the modulated state. Impor- 
tantly, the transition is first order on both sides of this 
point, i.e. we find first order transitions both in the uni- 
form and into the LO state. At lower temperatures the 
transition to the inhomogeneous superconducting state 
is second order. This is natural within our picture since 
the thermal fluctuations "die out" as the temperature is 
lowered. Within the present framework we cannot de- 
termine whether, should the quantum dynamics of the 
magnetization be accounted for, the transition would re- 
main first order to the lowest temperatures. However, 
since t*^Q ~ OMc for d-wave order parameter (Figs. [2] 
and [3]), it appears likely that the LO transition becomes 
second order again at high enough temperatures so that 
the quantum fluctuations are unlikely to have a major ef- 
fect. We now investigate the thermodynamic signatures 
of these transitions. 
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FIG. 6: (Color online.) Optimal wave-vector at N-SC tran- 
sition for s-wave with f/x ~ 0.4. Small-g approximation pre- 
dicts tLo,grad and t^^Q .grad- Upswing in q below t*^o,grad in- 
dicates breakdown of small-g approximation^. 



B. Thermodynamics at N-SC transition 

1. Specific heat jump at the second order transition. 

The specific heat jump, measured relative to the BCS 
s-wave value, at the second order N-SC transition along 
Bc{T) is given by 



^C/TrM 7C(3) {a'f 
1.43C/r, 



cO 



87r2 



27 



(33) 



Here again the prime denotes the temperature derivative, 
and the quadratic and quartic coefficients are determined 
from Eqs. (|2T|) and (f24| evaluated at the transition point 
and optimal modulation vector Qq. The results are pre- 
sented in Fig. [7] for the s and d-wave superconductors. 

Not surprisingly, the specific heat jump diverges on 
approaching the first order transition region. Note that 
in the absence of fluctuations, even though the transition 
remains second order throughout, there is a singularity in 
AC/Tc due to the vanishing of the quartic coefficient at 
Tp. The shoulder in the specific heat in the modulated 
state is found both with and without coupling to the 
magnetic moment, and hence simply refiects the details 
of the variation of the coefficients and the modulation 
wave vector with temperature. 

Of more interest is the low temperature behavior. 
While for s-wave superconductors the specific heat jump 
vanishes as T ^ for both 77 = and 77 ^ 0, for the 
d-wave symmetry the same jump is a) finite for 77 7^ 0, 
and b) exhibits a minimum at the lowest T. 

The key to understanding this behavior is in evalu- 
ating the T = limit of the coefficients aLo and 7lOj 
which can be done analytically as detailed in Appendix [Cl 
Note that the classical fluctuations disappear at T = 0, 
as evidenced by the linear in t fluctuation corrections 
in Eqs. (PT|) and and that the values of be and Qq 




T /T 



cO 



FIG. 7: (Color online.) Specific heat jump at second order N- 
SC transition for (a) s-wave and (b) d-wave symmetries with 
fix =0.0, 0.3. Inset: Specific heat jump AC at second order 
N-SC transition. 



at t = do not depend on r/ or x- For s-wave sym- 
metry, the optimal wave vector and critical fleld are 



e~^^^o^ ~ O.SeCcT^ and 6c,s = e"'^ /2 



0.28 



(7^; « 0.577 is Euler's constant) at zero temperature, 
respectively. We flnd that for the s-wave case in the ab- 
sence of fluctuations at Qo,s, the quartic coefficient jlo 
diverges as (6^ — (Qo,s/2)^)^'^/^ as the field approaches 
bc.s (see Eq. (|C10[) ). Hence AC/Tc = at zero tempera- 
ture irrespective of the value of 77. 

In contrast, we find that at zero temperature the opti- 
mal wave vector for the d-wave gap is 



Qo,i 



e '"^ exp 



V3- 1 



^0-1 « 0.67Co-i 



(34) 



with ba4/{2QQ^dio) ^ (1 + V^)/4)i/2 Si 0.83, and the co- 
efficient ^Lo{T = 0) = 0.07 remains finite for all values of 
77. The vanishing of the specific heat jump in the absence 
of magnetic fluctuations is now due to the vanishing of 
a' at r = (discussed in Appendix IC 2\ . The tempera- 
ture slope of the quadratic term, q;'(0) = a'^Q(O) -|-77x(0), 
increases as 77 becomes flnite, and this leads to a finite 
value of AC/Tc for d wave order in the limit T — in 
the presence of the fluctuations. 

The negative slope ait — Qoi the speciflc heat jump for 
d-wave (Fig.[7Jb)), is due to the temperature dependence 
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FIG. 8: (Color online.) Thermodynamics at first order N-SC 
transition for s-wave with fjx =0.28, 0.32, and 0.36. Shown 
are (a) decrease in entropy and (b) decrease in susceptibility. 
Temperatures t*p, Ilo, and t|,o are labeled for fjx = 0.28. 



of X- To explain this, we expand Eq. (j33p in t to find 



AC/T, 7C(3) / {a'Y a' {2^a" - fa') 



AC/T, 



cO 



27 



272 



(35) 



where all the derivatives and 7 are evaluated at i = 0. 

As discussed in Appendix IC 21 for d-wave symmetry 
the quadratic derivative q;'(0) — ?yx(0) is positive while 
the quartic derivative 7'(0) = — (7yx(0))^/(6iF) is neg- 
ative at low t. The second order quadratic derivative 
is a"(0) = a'/o(0) + 277x'(0) with a'/o(0) « 4.54 and 
x' = — x/^s/ is always negative for x{T)- Hence, with 
x(r), the initial slope at f = of the specific heat 
jump is determined by how strongly the fluctuations com- 
pete with superconductivity. As shown for fjx — 0.3 in 
Fig.[71[b), moderate coupling is sufficient to make promi- 
nent the dip in the specific heat jump for d-wave at 
low temperatures. For constant susceptibility, however, 
x' = 0, and the specific heat jump always increases from 
its value at i = 0. 



2. Entropy and magnetization at the first order transition 



Between t*p and ^^q, where the transition is first order, 
we compute the entropy jump, AS' = —df/dt, at the 



transition, and show it in Fig. ^a). From Eq. ([9]), the 
entropy jump is 



-AS = 



da 



\So\ 



at 

da 7 



dt 2v 



dt \2iy) 



dv 
'dt 



du r ^ y 
'dt \2^) 



(36) 



with |(5op = —^l(2v) at the first order transition. As 
the effective coupling parameter between magnetism and 
superconductivity, 77X: grows, more and more entropy is 
transferred at Tc from the magnetic fluctuations to super- 
conductivity, and the entropy jump increases. We find, 
as expected, that AS* is largest in the vicinity of ^lOj 
where (5o takes its maximum value, and is on the order of 
a few percent of the entropy difference between the SC 
state at T = and the normal state at T^iB — 0). We 
also find that 8q is moderate at the first order transition, 
with its largest value 5o(^lo) ~ 0.3(5o(i = 0, 6 = 0), and 
the results of our small (5o expansion make physical sense. 

The mismatch in the entropy jump in Fig.[8fa) at tj^o 
results from averaging the LO gap amplitude over the 
system size in the limit g = 0. Near iLO? the wavelength 
Afflo — ''^S,i)/q of A(a;) becomes comparable to the sys- 
tem size, and the profile of the order parameter near 
t^o resembles a single kinbiiii^ profile that describes the 
uniform-modulated transition within the SC phase. Be- 
low t^o the modulation vector go rises rapidly along bc{t), 
and the spatial averaging of the order parameter is jus- 
tified away from the immediate vicinity of Ilo- There- 
fore we expect that a calculation free of the single-mode 
ansatz, will give a greater entropy jump in modulated 
state in the immediate vicinity of Tlo- 

Since at the first order transition 7 changes sign, we 
expand this coefficient near t*p and t*j^Q along the transi- 
tion fine, 7 = gi{tc — t*) where gi is positive (negative) 
near ti = t*p (t^Q)- We find that, near the tricritical 
points. 



AS^SNiQ-Sscitc) 



-9i da 
2v dt 



(tc-t*), (37) 



where Sn and Ssc are the entropy in the normal and 
SC states, respectively. Hence —AS increases linearly in 
tc — ti as seen in Fig. [Sl^a). This behavior may be tested 
experimentally in magnetocaloric measurements. 

Exactly at the points t*p and i^Q, the entropy differ- 
ence between the normal and the superconducting states 
is zero. Instead, there is a rapid release of entropy upon 
lowering the temperature at a fixed field, and entering the 
SC state. Near tp and i^^, both the quadratic and quar- 
tic Landau coefficients are small and can be expanded 
about t*, namely a = ai{t — t*) and 7 — gi{t — t*). We 
then find 



\5.? = 



-g,{t ~ t*) ± ^gl{t - tXf - ia^vW^) 



(38) 



and, sufficiently close to t^, |(5o P is dominated by the tem- 
perature dependence of second term under the square 
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root. Thus, the entropy relative to the normal state 
varies with temperature as 

Sscit) ~ S^{t) ^ -^yS^tf-TI (39) 

where Sscit) and SN{t) are the entropy in the SC and 
N states, respectively. 

To further test the validity of our parameter choices, 
we calculate 

M/H={x-'+v\^o\Ley\ (40) 

along the first order N-SC transition. Here lAol^^e is 
the spatial average of the SC order parameter. We find 
that magnetization is suppressed by the onset of super- 
conducting order (see Eq. (jl5l) ) as entropy is transferred 
between their respective degrees of freedom. The frac- 
tional change in magnetization is 

AM _ TO(T)|Ao|i, 
M l + r7x(r)|Ao|Le ^ ' 

across the transition. This jump, as shown in Fig. [SJb), 
resembles the entropy jump in Fig. EJa), which makes 
sense as both quantities depend on the value |Ao| takes 
upon entering the SC state. Thus, the jump in So 
across the transition may be revealed by measuring both 
AM/M and AS* along the line of first order transition. 

Since we include fluctuations phenomenologically, it is 
possible that the first order transitions are due to an un- 
reasonable choice of the coupling parameter 77 such that 
the magnetization is strongly renormalized. As a check 
on the validity of our model, we verify that the magneti- 
zation does not change drastically at the N-SC transition. 
As shown in Fig. [SKb), the relative change in M/H at a 
first order transition is generally less than a few percent 
and validates our method of including of the magnetic 
fluctuations. 



While we considered only classical thermal fluctua- 
tions, since the regime of interest occurs for tempera- 
tures T/Tc ^ 0.2 — 0.5 we believe that this approach 
is sufficient. The question of whether the transition re- 
mains second order as T — > 0, e.g., when accounting for 
the quantum dynamics of spins, is left for future studies. 
Among other potentially interesting issues are whether 
impurity scattering, which is known to suppress the in- 
homogeneous LO state, enhances or shrinks the first or- 
der transition regime, what the results of combining the 
Zeeman field with the orbital coupling and vortex physics 
would be. In the present form our approach outlines a 
new, generic, path towards a first order N-SC transition, 
and demonstrates one experimentally observed feature: 
the separation between the onset of the first order tran- 
sition and the transition into a modulated state. It sug- 
gests that accounting for magnetic fluctuations which are 
known to exist in heavy fermion and other related com- 
pounds affects the shape of the transition lines, the order 
of the transition, and the behavior of the thermodynamic 
properties at the transition. 
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Appendix A: Expansion in Aq 

To derive the coefficients for the Landau free energy 
functional, we begin with the requirement = —Thi{Z) 
is an extremum with respect to Aq and A* so that 



V. CONCLUSIONS 

Motivated in part by experiments on the 115 heavy 
fermion compounds we considered the effect of a Pauli- 
limiting Zeeman field on s- and d- wave supercon- 
ductors in the presence of classical magnetic fiuctua- 
tions. We considered both the uniform and inhomoge- 
neous (FFLO) superconducting states, found that the 
amplitude-modulated state with the modulation vector 
along the gap nodes is favorable in the c?-wave case, and 
investigated the order of the transition. Our main finding 
is that there exists a range of temperatures, in the vicin- 
ity of the onset of the modulated state, where the normal 
metal-to-superconductor transition is first order both into 
a uniform and into a modulated state. The width of the 
temperature range increases with the strength of cou- 
pling to the magnetic fluctuation and is generally greater 
for s-wave systems. 



(Al) 

We construct the Landau free energy functional by ex- 
panding in powers of Aq, but not in the modulation wave 
vector q, which allows us to treat the low temperature 
region. To carry out this expansion we use the Gor'kov 
formulation of the Green's function approach. The nor- 
mal, 

G.(k,k';r) = - {Tr (ck,.(T)4, ,,(0)) ) , (A2) 
and anomalous, 

i^t(k,k';r) = - {Tr (i-(r)4,^+(0))) , (A3) 
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Green's function satisfy 

{iujn - ek,+)G'+(k, k'; iiu„) 

+ ^3;(k)AqFt(-k + q,k';i^„) = 5k,k', (A4) 
q 

(iw„ + e_k,-)^''''(-k,k';ia;„) 

+ ^3;(k)A;G+(k + q,k';iwO =0, (A5) 
q 

respectively. Here denotes imaginary time ordering, 
and ujn = 2itT (n + j) is the fermionic Matsubara fre- 
quency. The thermal average entering the free energy 
expression, Eq. (lAip . is given by 



where 



G'j_(k; iujn) = [ii^n - ek,±]' 



(All) 



is the normal state propagator for an electron in a Zee- 
man field. After integration over k, the interaction 
strength |A| in the quadratic coefficient a will be elim- 
inated in favor of the zero-field transition temperature 
TcQ. Assuming a circular Fermi surface, we use 2D angu- 
lar basis functions 

{1 s— wave 

V2 (kl - kl) = V2cos(20) d,2_y., 



-k,-4+q.+ ) = -T^i^^(-k,k + q;iL.„) . (A6) normalized so that {\y{9)\ 



/ FS 



We iteratively expand Eqs. (jA4p and (jA5l) in powers 
of Aq and A*, and hence find the series expansion for 
k, k'; icLi„). Using this expansion for the thermal 
average in Eq. (|A1|) . we integrate term by term with re- 
spect to Aq, and we obtain the Landau free energy den- 
sity, Fl — T jl? where I? is the 2D system size, up to 
0(|Ao|S) inclusive. We find 

F,. = ^5q,|Aqj2 

+ 7qi,...,q4^qi^q2^q3^q4^qi+q3,q2+q4 

{q.} (A7) 
+ <. -q6AqiA;^Aq3A;^Aq,A;^ 

{q.} 

^ '^qi+q3+q5,q2+q4+q6 

where the summation over {q^} includes all possible com- 
binations of the allowed Fourier components of A(r). 

The fully q-dependent coefficients of the Ginzburg- 
Landau expansion are given by 

"1 = m - |3^(k)pG"+(k + q;ic.„) 

„,k (A8) 

^0 



xG'i(-k;-ic^„), 



T 



7qi.....q4 = ^ E (k + qi; zc.„) 

n,k 

X G° (-k + q3 - q2;-«w„) 

X G°.(k-Hq2;iw„)G° (-k;-iw„), 



(A9) 



T 



E|3;(k)rG"+(k + qi;zc.„) 



X G" (-k + qs + q5 - q2 - q4; -»w„) 

X G° (k -f q2 + q4 - Qs; «w„) 

X G" (-k + qa -q2;-iw„) 

X G']_(k + q2;ia;„)G°_(-k;-ia;„), 

(AlO) 



1 where d is the az- 
imuthal angle in momentum space. Here (• • ■)pg indi- 
cates an average over the 2D Fermi surface. All the 
momentum sums are evaluated using the fact the the 
Green's functions are peaked at the Fermi energy, so that, 
for our model of a 2D circular Fermi surface 



k 



"2^ 



de. 



(A13) 



Appendix B: Tridiagonal Integral 

For the case of single-mode cos(Q • r) modulation of 
the order parameter, the magnetic contribution to the 
free energy functional (due to the off-diagonal k, k ± 2Q 
coupling) takes the tridiagonal form 

J-(M(r)) =TE(akMk-M* 

k (Bl) 
+ 6kMk-(M*+2Q + M*_2Q)) 

where 

ak = a^^Q + ir?|Ao|^), Vk (B2) 



and 



&k = &^-^^|Aop, Vk. 



(B3) 



This yields the partition sum 

k 

+ 6Mk-(M^+2Q + MUQ))] 
= nn / 2'(^^M)exp [- (a|Mk,.p 

1=1 k 

+ 6Afk,.(Af^+2Q_, + Af*_2Q,,))] 



(B4) 



where the product over i accounts for the d spatial com- 
ponents of M(r). To compute this integral, we separate 
the product over all wave vectors into a product over 
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components parallel and perpendicular to the direction 
of Q. As the terms comprising Z have no functional 
dependence on i, we have 2, 



Zq where 



k± k 

X exp ^bMk^,k 



{ml 



(B5) 



,fc||+2Q + ^^fcj_,/C|,-2g 



Due to the coupling between Mk^^k^^ and Mk^^kii±2Q, 
the product over fc|| can be divided up into a product of 
integrals taken only over wave vectors \k^^\ < \Q\, effec- 
tively employing the Brillouin zone method of solid state 
physics with |A;||| < \Q\ corresponding to the first Bril- 
louin zone. Each term in the product over < \Q\ is 
then an integral connecting fc|| to fc„ = + n(2Q), where 
n is an integer. The sum over k is cut off at a wave vector 
on the order of the inverse lattice spacing kc = tt/I. So, 
to cut off the sum over n, we define the cut off integer Uc 
such that k±n^ = fc|| ± nc(2Q) « kc- 

Separating the product over fc|| in this way, and intro- 
ducing the notational shorthand 

M„(fc_L,fc||) = Mfe^_fe||+„(2Q) , 

our partition sum can be rewritten 



^0 = n n y [• • • lyiMiMM^iMMo)] 

X exp [-a|MoP - bMo (M^ + M*_^) 



(B6) 



We integrate over the real and imaginary parts of M„ = 
M'^ + iM^ and restrict the product over k to be over 
one-half of fc-space because Mk — M^^^ for real M(r). 
However, as the integrand factors into two identical inte- 
grals over and M", we can integrate over AI^ alone 
and take the product over all values of k. Thus 



^o-n n 



■V{M[)V{M'_i)V{M[,)] 



(B7) 



X exp [-a(Mo)^ - 26Mo {M[ + MLi) - 



Beginning with n — 0, we integrate recursively over 
all and denote by a„ and the renormalized coef- 
ficients of (Af4n)^ and M'^M'_^, respectively. Working 
with a and b given in Eqs. (|B2I) and (jB3p . the integration 
coefficients are 



ao — a and bo — b for n — 0, 



^2 ^2 

fli = oq and bi ~ — for n — 1, 



and, for n > 1, the remaining terms 
a„6o 



''n+l 



bl 



and bn+i 



62 



(B8) 



(B9) 



(BIO) 



are determined recursively. The partition sum becomes 
Zo = 



A; I .k 




(Bll) 



k\,k 



where fc|[ e (— Q, Q) is understood. The free energy func- 
tional J^(M(r)) can now be replaced with its thermody- 
namic average J^ = — /3~^ln(Z) which is 



d_ 

2P 



|k|=0 



k I ,fc|i 



, (B12) 



where for the second sum n G (— nc,nc)- The prime 
implies that n = is excluded from the sum since the 
n = term is ln(l) = 0. 

In order to obtain the necessary small A expansion of 
Eq. (|BT2)) . we need to expand (a^ - bl) /al to 0(|Ao|^) 
inclusive. We do this by introducing recursion relations 



ao 



'1, n>l 



(B13) 



dn = an - bn = oq - "T" , n> I 



with the initial values si = oq and di = qq — 26g/ao, 
respectively. Taking qq — a and bo = b from Eqs. (|B2p 
and (jB3p . we expand Sndn/aQ = (a^ — b^) /oq to third 
order in b since b oc jAgp. Expressing a and b in units of 
l/2xT, we have the initial values 



si = a = I — 2b 
262 

di — a 



= 1 - 26 - 262-46^-^0(6'^), 



1 - 26 

and the remaining terms for n > 1 

s„ = (i„ = 1 - 26 - 62 - 26^ + 0(6'') 



(B14) 



(B15) 



With these expressions, wc find that s„(i„ — 1 — 46 + 262 
is independent of n when expanded to third order in 6. 
Thus, 



Sndn 1 - 46 + 262 



= 1 - 262 - 



+ 0(6^), (B16) 



al (1-26)2 
and, substituting 6 = — r7x|Ao|2/4, we finally obtain 

= 1 - lv\Mo\' + Wx'\Ao\'. (B17) 



up to 0(|Ao|^) inclusive. Since the summands no longer 
depends on n, we recollect the summation over k±, 
and n into a sum over |k| < k^- We take the sum to 
include all n e (—nc,nc) with the n = term identical 
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to the rest. We justify this by noting that, for a system 
of size L^, the sum over fc|| for n = is of order 2QL and 
is much smaller than the sum over all < kc (of order 
2kcL) since Q <^ kc (where Q ^ £,q^ and kc = tt/^). 

After subtracting the average magnetic contribution 
to the normal state, the fluctuation contribution to the 
superconducting free energy is 



J'LO,M = 



2(3 



E 

|k|=0 



ln( l + i,7X|AoP 



(B18) 



ln(l-i^V|Aor + kx'|Ao|' 



which, with d = 3, is the expression given in Eq. ([23 



Appendix C: Zero temperature limit 

We determine the Landau coefficients a^o ^-nd jlq 
and their temperature derivatives in the limit i — >■ 0. We 
first derive analytically aLO and jlo from Eqs. ([6]) (with 
prefactor of 1/2 for LO) and ([7]) for both s- and d-wave 
at zero temperature; there the magnetic fiuctuations die 
out so that ulo — c^lo and jlo = ILO- We then de- 
termine numerically the derivatives Q!'^q(0), Q!'/q(0), and 
7lo(0) for d-wave symmetry and add to them the mag- 
netic fluctuation corrections at i = 0. 



1. Analytic determination of be, qo, and 7lo 

We first evaluate the quadratic Landau coefficient for 
the LO gap modulation. In the limit t ~ 0, the quadratic 
coefficient becomes 

aLO^l(\y{e)\'\n[ib + qf))-^^(^^^ , (CI) 

where y{e) = 1 and y{9) = ^2 cos26' for s- and d-wave 
gaps, respectively. Here we use b = fiB/{2'KTco) and 
q = qcos{9 — 6q) where q = £^oQ/2. The angle 9q is 
the modulation direction with respect to the crystalline 
a axis, and 6q = tt/A for nodally-oriented d-wave. Inte- 
gration over 9 yields 



1 



In 



2 I 2 



(C2) 



where 7^ 0.577 is Euler's constant, and for d-wave 



■ exp ( -20-^ + 2a2 - i ) ~ 0.337 , 



(C5) 



bc,d = aqo,d ^ 0.278 , 



where a = {{1 + y/i)/Af''^ ~ 0.826. 

To determine the quartic Landau coefficient, first 
note that, in the limit T = 0, the Matsubara sum 
27rr^^F(aj„) becomes the integral / dujF{uj). Thus, 
we rewrite Eq. ([7]) as 

7.0= Re/ -\yiO)\%i^A,,0) (C6) 

where lj = uj/{2t:Tco) and 

(cj + ib) (3 (w + ibf - 



I^{uj,b,q,9) 



((w + Ibf 



(C7) 



We perform the angular integration changing variables 
to z = e'^, and then integrating around the unit circle 
in the complex z-plane. After thus averaging over the 
Fermi surface, we arrive at 

2q^ + 5q^iDl + 6Qt 
I,A^,b,<l)- 2c^g(g^+c^g)5/2 (C8) 



and 



^ ,_ , , 24Cji,8q^Ldb + Uq^ui^ + 40Ldl 



2 I ,7,2 



^{q^ + 24q'Cjl+40Cot) 



(C9) 



for s- and d-wave respectively. Here ioi, = uj + ib. Evalu- 
ating the integral over ui we arrive at 



ILO,. 



3 36^ - 2g2 

32^5(62-52)3/2 ' 



(CIO) 



and 



1 / 52 / 52 54^ 

7LO,.-:,TZ777(l-2-(3-36-+4G-j) ) (CU) 



dAn^q"^ \ q^ 



the quartic Landau coefficients at t = 0. From Eq. (|C10|) . 
we find that 7lo,s diverges as be — > qo (see Eq. (|C4P ) 
while, from Eq. JCl]) and (jClip we see that 7^0, d - 0.070 
remains finite when T — >■ 0. 



and 

for s- and d-wave respectively. We locate the transition 
by finding the maximum b for which aLO,d = 0, and we 
find that for s-wave 



2e*(i/2) = ij^ ^ 0.281 



(C4) 



2. Evaluation of derivatives for d-wave at T = 

The temperature derivatives of the quadratic and quar- 
tic coefficients are 

aLo(^,i) = a'Lo{t) + \vX (C12a) 
«Lo('7,i) = <oit)+Vx' (C12b) 
l'Lo{v,t) = 7ioW-^^V (C12c) 
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FIG. 9: (Color online.) Temperature derivatives of quadratic 
and quartic Landau coefficients for d-wave at fixed bc,d and 
qo,d in zero temperature limit. Main figure shows q^q and 
~7l05 both of which limit to zero at t = 0. Inset: a'[o ~ 4.54 
at t = 0. 



Expressions for a'j^Q^t), a'lQ{t), and 7^0 (0 are obtained 
by taking the first and second derivatives of Eq. ^ (with 
prefactor of 1/2 for LO modulation) and the first deriva- 
tive of Eq. (O, respectively, with respect to t. 



We determine a'^Q(O), q;'/q(0), and 7^0 (0) numeri- 
cally by fixing b — hc^d and q =^ qo^d and evaluating the 
derivatives as t approaches zero. As shown in Fig. |9l 
we find that a' 



LO 



i'lo 



and a'^Q 



4.544 at 



t = 0. Using these values and working with x(^)> we 
obtain a'^o(O) = 6.18?7x, a'loi^) = 4.544 - 152.8^x, and 
7lo(0) — ~l-48(77x)^ in the zero temperature limit. 
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